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Abstract 

We start from the QED Lagrangian to describe a charged many-particle system cou- 
pled to the radiation field. A covariant density matrix approach to kinetic theory 
of QED plasmas, subjected to a strong external electro-magnetic field has recently 
been developed [1,2]. We use the hyperplane formalism in order to perform a man- 
ifest covariant quantization and to implement initial correlations to the solution of 
the Liouville-von Neumann equation. A perturbative expansion in orders of the fine 
structure constant for the correlation functions as well as the statistical operator is 
applied. The non-equilibrium state of the system is given within generalized linear 
response theory. Expressions for the susceptibility tensor, describing the plasma re- 
sponse, are calculated within different approximations, like the RPA approximation 
or considering collisions within the Born-approximation. In particular, the process 
of relativistic inverse bremsstrahlung in a plasma is discussed. 

Key words: relativistic kinetic theory; QED plasma; hyperplane formalism; 
inverse bremsstrahlung, relativistic linear response theory 



1 Introduction 



In recent years the theoretical study of dense relativistic plasmas is of in- 
creasing interest. Such plasmas are not only limited to astrophysics, but can 
nowadays be produced by high- intense short-pulse lasers [3-5]. In view of the 
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inertial confinement fusion, one has to consider a plasma under extreme con- 
ditions which is created by a strong external field. This new experimental 
progress needs a systematic approach based on quantum electrodynamics and 
methods of non-equilibrium statistical mechanics. 

Considerable attention has been focused on a kinetic approach, formulated 
for the fermionic Wigner function. Using the Wigner operator defined in four- 
dimensional momentum space [6-8], a manifestly covariant mean-field kinetic 
equation can be derived from the Heisenberg equations of motion for the field 
operators. In this approach, however, it is difficult to formulate an initial value 
problem for the kinetic equation since the four-dimensional Fourier transfor- 
mation of the covariant Wigner function includes integration of two-point 
correlation functions over time. Describing modern pump-and-probe laser ex- 
periments by an initial value problem can lead to a significant simplification, 
since the pump and probe process can be described separately. However, at 
short time scales the plasma response depends on its initial correlation, which 
have to be included appropriately. 

An initial value problem can be formulated within a one-time formulation, 
where the field operators are taken at the same time and only the spatial 
Fourier transformation is performed. In the context of QED, the one-time for- 
mulation was proposed by Bialynicki-Birula et al. [9] and used successfully in 
their study of the QED vacuum. Within this approach the one-time Wigner 
function has a direct physical interpretation and allows to calculate local ob- 
servables, such as the charge density and the current density. The description 
in terms of one-time quantities is quite natural in kinetic theory based on the 
von Neumann equation for the statistical operator and provides a consistent 
account of causality in collision integrals. 

It should be noted, however, that the one-time description does not contain 
the complete information about one-particle dynamics. Spectral properties of 
correlation functions are naturally described in terms of two-point Green's 
functions which are closely related to the covariant Wigner function. In non- 
relativistic kinetic theory, where two-time correlation functions can, in prin- 
ciple, be reconstructed from the one-time Wigner function by solving integral 
equations which follow from the Dyson equation for non-equilibrium Green's 
functions [10], this problem can be overcome. In a relativistic theory it was 
suggested [11] that an expansion in terms of energy- moments can recover the 
complete spectral information within a one-time formulation. The discussion 
of spectral properties will not be issued in this work. Recently the aspect 
of relativistic kinetic theory was studied within the mean-field approxima- 
tion [11,12]. 

In this work we follow the studies presented in [1,2] where a density matrix 
approach to kinetic theory of QED plasmas subjected to a strong electromag- 
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netic field was considered. The BGR scheme [9] was generalized in two as- 
pects. First, the one-time formalism was presented in a covariant form. This 
removes a drawback of the BGR theory which is not manifestly covariant. 
The covariant formulation will be performed using the hyperplane formalism. 
Second, we demonstrate how the non-equilibrium statistical operator can be 
expressed within the generalized linear response theory in a hyperplane for- 
malism. The linear response theory was successfully applied in many fields of 
modern physics and can serve as tool to investigate the rather complicated 
structure of general kinetic equations. Here we are interested in the response 
of the plasma to an external perturbation, caused by the laser pulse. Consid- 
ering the external field and the response of the system as a small perturbation, 
the non-equilibrium statistical operator is expanded up to linear order with 
respect to the external field. It is clear, that this approximation is only reli- 
able for moderate laser intensities, but will break down for strong fields. The 
method applied here was successfully used in non-relativistic theory [13,14]. 

Having determined the non-equilibrium statistical operator, response func- 
tions like the susceptibility tensor is given. As an illustration, the relativistic 
susceptibility tensor is derived in the mean-field approximation. In the present 
work we also give first results beyond the RPA where interaction processes be- 
tween the charged particles and the EM field are considered. The susceptibility 
tensor, given by the current-current correlation function, is most appropriately 
expressed in terms of the force-force correlation function, which allows for a 
well defined perturbation expansion with respect to the interaction. It is shown 
how the absorption coefficient is related to the imaginary part of the force-force 
correlation function in lowest order perturbation theory. In this approxima- 
tion the absorption of an external field is treated. As well known, this first 
order process is possible only in combination with the Coulomb interaction 
in order to obey conservation laws. In Born approximation, considered here, 
the inverse bremsstrahlung is obtained. We give explicit expressions for plas- 
mas near the equilibrium. A final expression for the inverse bremsstrahlung 
can be compared with other results, like relativistic S-matrix calculations in 
vacuum [15] or with non-relativistic results [16]. 

The paper is organized as follows. In Section 2 we demonstrate how the quan- 
tum Hamiltonian defined on a hyperplane can be obtained from the classical 
QED Lagrangian. Major points like the gauge-fixing or the canonical quan- 
tization on the hyperplane is reviewed. Further the relativistic Liouville-von 
Neumann equation is formulated on the plane. In Section 3 it is shown how the 
Liouville-von Neumann equation can be solved in linear response, fixing some 
initial distribution. The self-consistency relation leads to a response equation, 
which defines the four dimensional susceptibility tensor. The result for the sus- 
ceptibility tensor in RPA is presented. In section 4 the inclusion of collisions 
is discussed. The electron-ion collisions are considered in lowest order per- 
turbation theory. In this approximation the absorption coefficient for inverse 
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bremsstrahlung is related to the imaginary part of the force-force correlation 
function. Finally Section 5 concludes with a discussion of the results and gives 
a short outlook. 



2 Canonical Description in the Hyperplane Formalism 

In this section we demonstrate the derivation of a quantum Hamiltonian start- 
ing from the classical QED Lagrangian. The main issues like gauge fixing, or 
canonical quantization on the hyperplane are discussed to some extend. 

2. 1 The Lagrangian Formulation of the Plasma 

We consider a charge neutral system consisting of two species of fermions, 
like for instance electrons (mass m, charge Z e e, Z e = —1) and ions (mass rrii, 
charge Zie). In particular, we consider protons (Zi = 1), but the generaliza- 
tion to arbitrary charged particles is straight forward. The standard covariant 
formulation of a QED plasma is given in terms of the Lagrangian C 

C'(x) = C D (x) + C EM (x) + C mt (x) (2.1) 

where C D (x) is the Dirac part, describing the fermionic components of the 
plasma (e.g. electrons and protons), C EM (x) is the electro-magnetic part and 
C int (x) describes the interaction in the plasma. In standard notation [17,18] 
we can write the different terms as 

C D (x) = Y, f(x) <^rd, - m) r(x) , (2.2) 

C EM (x) = --^ v {x)F^{x) , (2.3) 
Ant(z) = " E Z c ej;(x)A»(x) , (2.4) 

c=e,i 

*r- > > < 

where cL = <9 M — <9 M - In the following the electro-magnetic field tensor is 
taken in the form = d^A u — dyA^. The current density four- vector will be 
expressed as j c = - ^ c 7 At V ;C with e < 0. 

Additionally to the terms given in Eq. (2.5) we consider the influence of an 
external field A^ xt (x) on the system. This external field is not a dynamical 
field variable, but is some given function of space and time. This implies that 
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there is no back reaction mechanism of the system onto A^ xt (x). The external 
field couples to the fermion current and therefore the complete Lagrangian 
takes the form 



C(x) = C'(x) + C cxt (x), 

AsctW = - 2 Z c e f»( X ) A ext(x)- 



(2.5) 
(2.6) 



c=e,i 



2.2 Gauge fixing 

The Lagrangian given in Eq. (2.5) is covariant and gauge invariant. However, 
as well known, the gauge invariance leads to non-physical degrees of freedom, 
which have to be eliminated from the description. One possible way to perform 
this elimination is to apply an additional gauge condition. We will follow this 
way and use the Coulomb gauge V • A = 0, which is the most natural choice 
for Coulomb systems, since the Coulomb interaction appears naturally in this 
description. 

The disadvantage of this method, however, is that the explicit covariance is 
lost. In particular, after application of some further approximations it is hard 
to control, if the final result will be expressible in covariant form. For this 
reason it is most advantageous to use the hyperplane formalism [19], which 
singles out time-like and space-like parts in a covariant way. 

Beside the natural appearance of the instantaneous Coulomb interaction due 
to the breaking of the explicit covariance by the Coulomb gauge constraint, 
there is a second advantage for this method. In view of laser plasma inter- 
actions in general we have to deal with a highly correlated plasma, which is 
subject to the external laser field. In order to describe the correlated plasma 
initially, the time coordinate has to be singled out, which leads to the break- 
ing of covariance. As already mentioned above, in the hyperplane formalism 
the inclusion of initial correlation as well as the gauge fixing leading to the 
Coulomb interaction can be performed in a manifest covariant manner. The 
main issues of this formalism will be discussed in the next section. 

2. 3 Introduction of space-like hyperplanes 

A space-like hyperplane a = o nr in Minkowski space can be characterized 
by a unit time-like normal vector n M and a scalar parameter r which may be 
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interpreted as an " invariant time" . The equation of the hyperplane cr reads 

x • n = t , n 2 = n^n^ = 1 . (2.7) 

In the special Lorentz frame where n M = (1, 0, 0, 0) and consequently Eq. (2.7) 
reads x° = r the parameter r coincides with the time variable t — x°. We will 
refer to this special frame as the "instant frame" , since only here observables 
are measured at the same instant of time t. Expressing the field variables 
as functionals of the hyperplane ^[(7 nr ] or as functions of n and r, The La- 
grangian can be expressed on the hyperplane u nT . The gauge condition is 
expressed on the plane according to 

= , (2.8) 

where the following decomposition of four vectors is used 

V = n»V ]} + V£, V }] = n v V, V? = A^V , (2.9) 

d , = + V„ V, = A;d u = A; A , (2.10) 

with the transverse projector A^ v 

A", = ^-n\. (2.11) 

Applying the decomposition (2.9) - (2.10) to the Lagrangian (2.5), we obtain 
* = lW - \ (VM, - A») (V^„ - A J) 

c 

where we have introduced the notation 

= WA" ± -V V A1 . (2.13) 

Due to the gauge condition (2.8) we find a constraint equation, similar to the 
Poisson equation 

WpE^'f. (2-14) 
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The solution of Eq. (2.14) is 

A ll (r,x ± )=Y / Z c e f da'G(x ± -x' ± )j;(r,x' ± ) , (2.15) 

c L 

where the Green function G(x ± ) satisfies the equation 

V^G(x ± )=5 3 (xJ (2.16) 

and the three-dimensional delta function on a hyperplane a n defined as 

5\ Xa ) = J*V- e -**6(p-n). (2.17) 

The solution of Eq. (2.16) for G(x ± ) is given by 



(2tt)3 



G ( x .) = - / 7^r, e~ ip - x S(p ■ n) - . (2.18) 



The variable A, can now be eliminated from the Lagrangian density (2.12) 
using Eq. (2.15) imposing appropriate boundary conditions. Then a straight- 
forward algebra leads to 



+ C D + C cxt -yXJ Z c Z c J da'j^r, xJG(x ± - <)jf(r, x'J .(2.19) 

C,C (j 



2.4 The commutation and anti- commutation relations 



In the canonical quantization scheme, which will be applied here, the commu- 
tation and anti-commutation relations of the canonical field operators have to 
be derived. We follow the method, originated by Dirac [20,21]. The momentum 
n ± canonical to the electro-magnetic field variable A ± ^ is defined as 



dC 

IT =^ = -A ln , (2.20) 
dA 
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with the r-derivative A ± . Similarly, for the fermionic field variables we define 
the canonical momenta 7f and it according to 

with the decomposition of Dirac's 7-matrices 

7 m = ^ 7|1 H+7^H, 7|l (n)=n^, 77» = ~ n»n v ) 7 ". (2.22) 

The dynamical fields and its canonical momenta will now be interpreted as 
operators, satisfying commutator and anti-commutator relations. These rela- 
tions can be obtained by calculating the Dirac brackets, which account for con- 
straints, like the gauge-fixing constraint Eq. (2.8). In Appendix A we shortly 
review the calculation, leading to the following non-vanishing relations 



A»(t,x ± ),W ± (t,x' ± )\ = tc^(x ± - x' ± ) , 



where 



( x± - x' ± ) = J 



d 4 p 
J2nf 



-ip-(x-x') 



S(p ■ n) 



l 1 V " 

p±p ± 
pi 



(2.23) 
(2.24) 

(2.25) 



For the Dirac field operators the anti-commutation relations on the hyperplane 
take the form 



^ac(^ x ±)^a'A T ^±)j = U\\( n )\ aa , 6 c,c>$ 3 ( x ± ~ x '±), ( 2 -26) 
^acC^J, $a'A T > X '±)\ = \iac( T , X J^a'c'( T , X 'j \ = 0, (2.27) 



where a, a' are the spinor indices and c, d denote the different species. In the 
special Lorentz frame where x^ = (t, r) and n M = (1, 0, 0, 0), we have 7,, = 7 
and 5 3 (x ± — x' ± ) = 5(r — r'), so that Eq. (2.26) reduces to the well-known 
anticommutation relation for the quantized Dirac field. 
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2.5 The Hamiltonian 



The quantum Hamiltonian can be constructed by a Legendre transformation, 
defined on the hyperplane as 

H(n)= J rfa|n ±M i^ + 7T^ + ^7r-£| , (2.28) 



where C is given by Eq. (2.19). The Hamiltonian can be written in the form 
= H D (n) + H EM (n) + H- mt (n) + ifj xt (n) , (2.29) 

where H D (n) and H EM (n) are the Hamiltonians for free fermions and the 
polarization EM field respectively, H- mt (n) is the interaction term, and H^ xt (n) 
describes the external EM field effects. In the Schrodinger picture the explicit 
expressions for these terms are 

H D {n) = E/ da$ e (-^7» V M +m c ) 4 , (2.30) 
H EM (n) = Jda (jF^F? - ^ft?) , (2.31) 
H- mt (n) = e —Y. Z c Z d J da j da'%(x ± )G(x ± -x' ± )jl (x' ± ) 

u n u n 

+ E [daZ c ejlAl, (2.32) 

c X 

n 

(n) = J2j daZ c e];(x ± )AUr,x ± ) . (2.33) 



#ext 



The field strength tensor F ± and the transverse field operators A± and ITf 
in Eq. (2.31) are defined according to their classical relations (2.9) and (2.13). 
The longitudinal part A,, has been eliminated in the interaction Hamilto- 
nian (2.32) by the operator version of Eq. (2.14). In Eqs. (2.30) - (2.33) normal 
ordering in operators is implied. The self-energy contribution to the last term 
in Eq. (2.32) is omitted, so that the product (x ± ): (x' ± ): is understood. 
The generalization of the Hamiltonian to a many-component case is obvious. 

It should be mentioned, that the same expressions for the Hamiltonian Eq.(2.30) 
- (2.33) can be obtained from the symmetrized energy-momentum tensor, the 
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so called Belinfante tensor T [22] via 

H(n) = PX = J danT^n" . (2.34) 



2.6 The relativistic von Neumann equation 



As already noted, the state of the system |\P[<7 nT ]) is taken as a functional 
of cr nr . A specific frame of reference can now be related to a family of space- 
like hyperplanes with a fixed normal vector n. The relation between different 
frames of reference is given by a homogeneous Lorentz transformation A 

°n,r <V, T = A(T n,r ■ X ^ x' = Ax . (2.35) 



Eq. (2.35) means that x^ = x'^n'^ = r, i.e. x is located at a and x' at a'. 

With a unitary representation of the homogeneous Lorentz group U (A) state 
vectors on different planes are related by [23] 

17(A) |tf[A<r]) = |*[<t]) . (2.36) 



Having fixed, the evolution of the state in this frame of reference is governed 
by a representation U(a) of time-like translations in the direction of the 
normal n^. The generator of this transformation is P M , the energy-momentum 
vector. U (a) is given by 

U(a) = exp{iP^} . (2.37) 



Writing the states as functions of n and r, we find for an infinitesimal time-like 
translation = 5t 

|*(n, r + St)) + i5r (P^) |*(n, r)> = |*(n, r)) , (2.38) 



from which we obtain the relativistic Schrodinger equation 

i^-\*(n,T)) = H(n)\*(n,T)) (2.39) 



with the Hamiltonian on the hyperplane given by 

H(n) = . (2.40) 
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In the presence of a prescribed external field, the energy-momentum vector 
and, consequently, the Hamiltonian H T (n) can depend explicitly on r. Com- 
bining Eq. (2.39) with the adjoint equation for the bra-vector, one finds that 
the statistical operator g(n, r) for a mixed quantum ensemble obeys the equa- 
tion 



dgjn, t) 
dr 



g(n,r),H T (n) 



= 



(2.41) 



which is analogous to the non-relativistic von Neumann equation. 

In order to solve Eq. (2.41), some boundary condition have to be imposed on 
the statistical operator. The standard boundary condition in kinetic theory 
is Bogoliubov's boundary condition of weakening of initial correlations which 
implies the uncoupling of all correlation functions to one-particle density ma- 
trices in the distant past, i.e., for r — > — oo. In the scheme developed by 
Zubarev (see, e.g., [14]), such boundary conditions can be included by using 
instead of Eq. (2.41) the equation with an infinitesimally small source term 



dg(n,T) _ . 
Or 



g(n,r),H T (n) = -r) {g(n, t) - g iel (n, r)} , (2.42) 



where 77-^+0 after the calculation of averages. Here g re \(n, r) is the so-called 
relevant statistical operator which describes a Gibbs state for some given non- 
equilibrium state variables. In QED kinetics these variables are the Wigner 
function and the photon density matrix. In general we will call these non- 
equilibrium state variables the relevant operators B^. In the next section we 
will search for solutions of Eq. (2.42) up to terms linear in B%. 



3 Linear Response Theory 



Using the results of the last section, it is possible to construct a kinetic theory 
on the hyperplanes, which was extensively discussed in [1]. 

Here we show how the non-equilibrium statistical operator for the case of small 
deviations from the equilibrium distribution can be constructed. This can be 
done by solving Eq. (2.42) with a appropriate choice for the relevant distribu- 
tion. This treatment, known as generalized linear response, was successfully 
applied in different studies [13]. 

Further for simplicity we consider the plasma in the adiabatic approximation, 
where the dynamics of the positively charged component of the plasma is 
frozen. 
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3.1 Fluctuations near equilibrium 



For a given relevant distribution g rel , a formal solution of the Zubarev equa- 
tion (2.42) is given by 

T 

g(n,r) = rj J dr' e" e ^ U(r,r') g iei (n,r') U\r,r') , (3.1) 



where the evolution operator can be written as the ordered exponent 

U(t, t') = T T exp { -i f H f (n) alf } . (3.2) 



After partial integration, the expression (3.1) becomes 

e (n, r) = g rel (n, r) + Ag(n, r) , 



(3.3) 



T 

Ag(n,r) = - J dr' e^"^ 



x U(t, t') 



dr' 



— i 



g iel (n,T'),H^(n)]\uHr,T') . (3.4) 



The Hamiltonian H, as defined in the Eqs. (2.29) - (2.33), is rewritten in the 
adiabatic approximation in decomposed form according to 



H(n) = H s (n) + H ext (n) , 

H s {n) = H D (n) + H EM (n) + H rad (n) + H int (n) . 



(3.5) 
(3.6) 



In (3.5) - (3.6) the system part of the Hamiltonian H s is decomposed into the 
kinetic parts H D and H EM and some part describing the collisions H int and 
H rad , which will be treated within perturbation theory. The different parts of 
the Hamiltonian are taken from Eqs. (2.29) - (2.33) 



H D (n) = I dai/j e [ ~H{n) V M +m e ) ^ e , 



(3.7) 
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H EM (n) = J da Qi^iT - ^ ± ^) , (3.8) 

H rad (n) = -eJda ? ±ll (x ± )AlM > (3-9) 

H int (n) = -efda (^(xJA^xJ +j?(xJA™(x x j) , (3.10) 

^cxtH = -ej da JiJCt^iJ . (3.11) 

In the following we will drop the index "e" and the electron spinors are denoted 
by ijj and ip. 

The radiation term, Eq. (3.9), is important for the description of photon emis- 
sion from the plasma or for photon scattering in the plasma. As an application 
we will focus to the absorption of a classical electro-magnetic wave by a rela- 
tivistic plasma in the next section, where the radiation term can be neglected. 

In what follows, the non-equilibrium state of the system must be specified. 
We use Zubarev's method [14] of a non-equilibrium statistical ensemble in 
linear response. The so called relevant statistical operator f? rcl (n, r) describes a 
generalized Gibbs distribution, which characterizes the initial non-equilibrium 
state of our system. In the linear response regime we will only consider small 
fluctuations from the equilibrium. The relevant statistical operator can be 
written in the form 

g iei (n, r) = Z~l((3, u, 0* ; r) exp ( - (3 itP^ - vQ 

(3.12) 



J da £, 



>1(x ± ;t)BZ(x a 



where the relevant observables B^(x ± ) define the non-equilibrium state and 
will be treated as small quantities. It should be mentioned, that the indepen- 
dence is carried completely by the set of Lagrange multipliers 0^(a; ± ; r). The 
first two terms in the exponential of Eq. (3.12) describe the generalization of 
the Gibbsian distribution g (n) of the grand canonical ensemble 



g (n) = Z-\P, v- n) exp j - (3 [n»P^ - uQ] } , 



(3.13) 



where 
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n"P„ = ff., Q = n t f . (3.14) 

The equilibrium statistical operator does only depend on the additive integrals 
of motion P , Q. 

In the following we demonstrate how correlation functions can be derived in 
linear response within the hyperplane formalism. This means, that all expres- 
sions will be approximated to first order in quantities, describing the deviation 
from the equilibrium. 

Expressing the relevant part of the statistical operator Eq. (3.12) in linear 
response, as demonstrated in Appendix B, we find in Fourier representation 



g vcl (n,r) = Q (n)+(3e- tk r T f dz 5>£(A;) B^(k ± ,izP)g (n) . (3.15) 

o * 

A similar calculation leads to the irrelevant part of the statistical operator in 
linear response as (see Appendix B for details) 

oo 1 

Ag(n,r)=(3e~ l ^ T J df e l(fc n +M?)f J dz 
o o 

£ [hY{-f + iz/3) - ik,BY{-f + iz/3)]<^(k) 

t 

+4L(k)l(-f + iz/3)}g (n) . (3.16) 

By definition, the mean values of the relevant observables are some prescribed 
functions of space and time and are obtained by averaging with the relevant 
statistical operator only. The self-consistency relations, by which the Lagrange 
multipliers are determined, can therefore be expressed as 

Tr {Bt(x ± )Ag(n, r)} = . (3.17) 

Multiplying Eq. (3.16) by B%,, taking the trace and rearranging the indices we 
find the response equation 



In Eq (3.18) the correlation functions (A;B) and its Laplace transforms 
(A;B) are defined according to 
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1 1 
(A-,b)= J dzTr{A(-if3z)B^g } = J dz Tr [AB\i(3z)q ) , (3.19) 
o o 

oo oo 

(i ; B) = J df e lT < f (i(f) ;B) = J df e ir < f (i ; B(-?j) . (3.20) 
o o 

Making use of Eqs. (3.15) and (3.16) local observables can be expressed in 
linear response, in particular the induced current j^ nd = 5(j^), is written as 

HO,) T = 5<0/» exp{-ifc„T}] 




(3.21) 



Eq. (3.21) describes the response of the system due to weak perturbations 
Ag Xt . The Lagrange multipliers can be eliminated from Eq. (3.21) by making 
use the self-consistency relation (3.18). 

If the induced current can be represented by a linear combination of the rele- 
vant operators = Y^v Q-i'B^ , Eq. (3.21) is simplified as 

where the remaining terms in Eq. (3.21) are canceled due to the self-consistency 
relation (3.18). The Lagrange multipliers can be eliminated in Eq. (3.22) by 
multiplying Eq. (3.18) by a e and summing over £. Further partial integration 
of correlation functions 

iz(A ; B) +(A;B)=(A\ b) = -(i ; B) (3.23) 
can be applied and finally we obtain 




(3.24) 



Now Eq. (3.24) can be plugged into Eq. (3.22) and we have the induced current 
expressed in terms of the external field within the linear response approxima- 
tion 
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= " + ih '^)} A ext • (3-26) 

In the second line again partial integration [Eq. (3.23)] was applied. 

The susceptibility tensor describing the response of the system to an 
external perturbation, is defined as 



3» d (k) = X^KxtW (3-27) 
and can be read from the equations (3.25) and (3.26) 

w = -P(i->L\ +lT1 = - iPKii ; ■ ( 3 -28) 

Eq. (3.28) gives the susceptibility tensor in terms of the current- current or the 
current-force correlation function. The major task to proceed is to evaluate 
these correlation functions within certain approximations. We will show how 
perturbation theory can be applied. The simplest evaluation of the correlation 
function is given in the RPA approximation, which will be considered in the 
next section. As a next step we demonstrate how collisions can be included. 

Since we are working in the adiabatic approximation, we chose a reference 
frame in which the the ions are at rest. This means that it is most convenient 
to use the instant frame formulation (n = (1,0,0,0)). 

3.2 The RPA-result for the correlation function 

The susceptibility tensor (3.28) can be most easily calculated in the RPA- 

approximation. In this approximation the force operator j u in Eq. (3.28) is 
calculated with the Dirac part in the Hamiltonian only. 

It is convenient to perform the calculation using the plane wave expansion for 
the field operators ip and ip according to 



= / jSkr 2 i^k ? t hpAPl s)e ~ iPX + 8)eHP ") ' (3 ' 29) 



(3.30) 
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We use the four-dimensional scalar product and the mass on-shell condition 



x ■ p = p ■ t — x ■ p 



Po = E P = +Jp 2 + m 2 . 



(3.31) 



The operators b, M, d and d) are the electron creation and annihilation oper- 
ators and the corresponding antiparticle creation and annihilation operators. 
The current operator can be expressed in terms of these operators 



r d 3 p pm . , J 

=e J (2^v^M 

/ -JH— :d ^ ,4 : (v(p + hk, s')7^(p, *)) e^'-W* 

V ^p+hk 

y 1 ^ -K-mAy- («(p - hk ' s 'hAP' s )) e-^-w* 



+ 



+4 



+ 4 



->p~hk 



m 



E 



p~hk 



m 



E 



p+hk 



■d- p+hk J P:S : (v(-p + hk, s') % u(p, sj) e-i'W 

A : " Hk > S 'K V & s )) e*<^ +J W*l. (3.32) 



The result (3.32) was obtained after integration over x and p'. 

In a similar way the Dirac Hamiltonian is written in the plane-wave represen- 
tation 



H D = J d 3 p £ E p (d\ s d ps + bl s b ps ) , 



(3.33) 



where all products of operators are taken in normal order. 



From the Heisenberg equation of motion and the current operator (3.32) the 
force operator is found 



re 



d 3 p 



m 



E 



m 



h J (2nhf \jE p ^ 

dlJ p+ m,s' (y{p + hk, s')^v(p, s)) 
e ^( E P - E p+hk) t 



E 



p+hk 



Ep - E p+hk 
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+ 



+ 



+ 



m 



E. 



p—Hk 



E p Ep-hk 



e -i( E P~ E P -hk) t 



m 



E, 



p—hk 



E P + E P -hk 



m 



E 



p+hk 



E p + E p+hk 



Q-Ti( E P+ E v -hk) t 

dl s bl p _ mjS , (u(-p - hk, s')7^(p, s)) 
e i( E p+ E v+ hk> t 



(3.34) 



Eq. (3.32) and (3.34) can now be used to calculate the correlation function 

uj+ir) ■ The calculation is shown in Appendix C. Together with 
the Eqs. (C.4) - (C.7) as well as Eq. (3.28) we have the 4-dimensional RPA 
susceptibility tensor in the form 



e 2 f 

x»A k ,u) = ~i J 



d 3 p 1 



(27rhr e p 



-1 


f(Ep) - 


f(E p+hk ) 


E p+ hk 


frw + E p - 


E p+hk + it] 


1 


f(E p ) - 


f{E p -nk) 


+ F 


hu-E p + E p _ m + if] 


1 


1 " f(E p ) 


~ f(Ep-Hk) 


Ep-hk 


hu-E p - 


■ E p _ m + 17] 


-1 


1 " f(E p ) 


- f(E p+hk ) 



■ E p+m ^ + E r + Ep+m + ir, - -M-X - ^ + •»]}} ^ 

The fermion and anti-fermion distribution functions, / and / respectively, are 
defined by Eq. (C.3). The equation (3.35) can be written in a more compact 
form, if we notice, that the susceptibility tensor can be decomposed as 



Xij 



k Q 



Xio 



~Xoi 



Xoi XiO 7 X i 

Kq 



k 1 



k 2 



1 3 



X 



, k { kj 

A ^2 Ajj 



X 



X 



ij ' 



(3.36) 
(3.37) 



This means that the longitudinal and transverse component completely deter- 
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mine the susceptibility tensor. This decomposition can be shown, making use 
of the current conservation uj — k ■ j = . 

In order to give an expression for the x l an d X* we calculate the Dirac traces 
of the spatial part x tj ■ For instance the first trace in Eq. (3.35) yields 



= 4{AGf + hf) j + (j + hf)^ 3 + 5 t] [E p E p+hk - p(p + hk) - m 2 ]} .(3.38) 

Calculating the remaining traces accordingly and projecting the longitudinal 
and transverse part [see Eq. (3.37)] as well as shifting p — > —p in the first and 
fourth term in Eq. (3.35), we find 



d 3 p j 

f(E p ) - f(E p _ hk ) 
hu-E p + E p _ hk + it] 



X i ' t (k,uj)=e 2 J 
A* 



i - Re v ) - m 



p—hk) 



frw + E p + E p _ hk + irj 



f(E p ) - f(E 



p—hk) 



hw + E p - E p _ m + ir]_ 



1 - f(E ) - f(E 



p—hk) 



hw-E„- E. 



p—hk 



+ if] 



(3.39) 



with the longitudinal and transverse projectors A 



El + hpk-2^- 

Ai = 1 ± — k — , (3.40) 

^p^p-hk 

m 2 -hpk+ { ^f 

A* ± = 1 ± — • (3-41) 

Observing, that A±* is invariant under the shift p — > —p + hk, we can perform 
this shift in all terms containing f(E p _ hk ) and f(E p _ hk ) in Eq. (3.39). We 
finally have the result for the longitudinal and transverse RPA susceptibility 
tensor 



d3 P 771/ \ f #p - E P -hk A e,t 



(2nh)* " p, \(Ep-Ep-hk) 2 -^-^V 

A?l+x2o (3-42) 



£p + E p _ hk , , , , , 



(E p + E p ^f -d*-iri 
F{E p )=2(f(E p ) + f{E p j) , (3.43) 
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( 



fc> W ) = _2e 2 / Ep + , A? • (3-44) 



It should be noted that the result (3.42) corresponds to familiar expressions, 
for example published by Tsytovich [24], given in the form of the dielectric 
tensor. 

The major difficulties in the theory arise, by the inclusion of collisions, which 
will be the main issue of the next section. 



3. 3 Inclusion of collisions 



As in the case of the RPA approximation, we restrict ourselves to the instant 
frame description, where the ions stay at rest (n^ = (1,0,0,0)). This also 
implies that we only have to consider the Coulomb interaction term of the 
electrons with the ions, which will be treated perturbatively. In the case of an 
arbitrary reference frame, where the ions are moving, we have n M = (n°, n) = 
(7, vy) , where 7 is the relativistic factor 7" 1 = — v 2 . 

Collisions are included into the response function Eq. (3.28) most conveniently, 
by making use of partial integration and writing Xliv 1X1 terms of correlation 

functions including the force operator j v in the denominator. We can find the 
matrix equation 



J Lb 1 Jv) 



(jn'Jp 

(j a 5 3u) M a/3 



M. 



a/3 



(3.45) 



M. 



(j? ; h) 

\3 a i Jv) \3 a i 3/3/ 



(j a 'iJ/3/ 



(3.46) 



In Eq. (3.46) the force-force correlation function appears, which is well suited 
for a perturbative treatment. Different processes can be found in each order of 
a perturbative expansion. In QED they are well studied and can be represented 
by Feynman diagrams. In the Figs. (1) and (2) we have shown the second 
order diagrams for the bremsstrahlung and inverse bremsstrahlung process. 
These processes are well known and expressions were given by QED S-matrix 



20 



■A 



Fig. 1. Feynman diagrams for second 
order bremsstrahlung of electrons in 
the Coulomb potential of the ions. 



Fig. 2. Feynman diagrams for second 
order inverse bremsstrahlung of elec- 
trons in the Coulomb potential of the 



ions. 



calculations [15] as well as by a non-relativistic treatment [16]. We will draw 
our attention to the inverse bremsstrahlung within our formalism in the next 
section. 



4 Inverse Bremsstrahlung 

As a second order process we consider in this section the absorption of electro- 
magnetic waves in a relativistic plasma. For illustration we derive the corre- 
sponding absorption coefficient in the long wavelength limit. 

4-1 The absorption coefficient 

Before calculating the inverse bremsstrahlung in second order, we need to give 
some relations, which allow to extract the absorption coefficient a out of the 
4-dimensional susceptibility tensor. 

In Eq. (3.36) and (3.37) a general decomposition of the susceptibility tensor 
into its longitudinal and transverse part was shown. The transverse compo- 
nent x* can a l so be projected out of by the two orthogonal transverse 
polarization vectors e (k, 1) or e (k, 2) 



While the susceptibility tensor describes the response of the plasma due to 
an external field, the dielectric tensor e connects the internal and external 
fields with each other. The longitudinal and transverse part of the dielectric 
tensor is related to the susceptibility according to [25] 



e^(k,i) X (k)^(k,i) 



i = 1,2 



(4.47) 




(4.48) 
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In the long wavelength limit, which will be taken here, it is seen from Eqs. (4.48) 
and (4.49) that the longitudinal and transverse part of the dielectric tensor 
and the susceptibility tensor coincide. The long wavelength approximation is 
well justified in the optical regime. The absorption coefficient is related to 
the damping of electro-magnetic waves, which is determined by the imaginary 
part of the dispersion equation. In the long wavelength limit we have 



a(u) = — — lim Im eHu, k) , (4.50) 
n(uj) fc^o 

where in the following the index of refraction is approximated by n(u) pa 1. 
From Eq. (4.49) we can find the relation 



lim ImeUu) = lim \ 2 lmX " 2 , (4.51) 

fc-^o v ; fc^o e 2 fc 4 r - 12 r - 12 '- 2 - 



— 4- 

o4 I 



Rex* + Imx 4 -2$-Rex* 



with a fc 4 -term dropped due to the limes k — > 0. It can be shown, that non- 
trivial solution for the transverse susceptibility, making use of Eq. (3.28) as 
well as Eqs. (3.45) and (3.46), can be written as 

X\k,u) = -P(j ± ;j ± ) v rT^-n '— ■ (4.52) 

-iuj (j ± ; j ± ) + (j ± ; 



In these expressions we use the transverse current j ± = e^j^. Contributions 

from the determinant, containing (j a ;j v ) are neglected, since these terms 
are of higher order. In order to derive Eq. (4.52), the current conservation, 
leading to the composition (3.36) and (3.37) was used. Treating the collisions, 

i.e. (j ± ; in Eq. (4.52) as small contributions, we can write to first order 
in the force-force correlation function 



(4.53) 



UJ 



Finally with Eq. (4.50), (4.51) and (4.53) we obtain the absorption coefficient 
containing collisions in first order 



aM = -/31im| I Re(i x ;j x )^. (4.54) 
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The force-force correlation function in Eq. (4.54) contains both, the emission 
and the absorption of photons. The evaluation of the correlation functions (see 
next section) shows that the absorption is described by 

<*(«) = -/? Jim ^Re(i ± (-k)-j ± (-k))_^ v , (4.55) 

with u > 0, k > 0. The emission of photons, which will not be considered here 
further, is related to the absorption by the interchange uj — > — uj and k — > —k. 

The relation, given in this section can straightforwardly generalized to arbi- 
trary hyperplanes. Since we work in the adiabatic approximation all expres- 
sions were given in the instant frame, which is off course most convenient 
here. 

4-2 Evaluation of the correlation function 

Evaluating the force-force correlation function in Eq. (4.55) we will follow the 
same way presented in the last section. The perturbative part in the Hamil- 
tonian, i.e. the Coulomb interaction with the ions, is written as 

H ion = ejd 3 x :j (x): A° lon (x) 

,/^:rf; i _„, r ,4 i , r :(0(p 1 -9,r')l 0l .(p 1 ,r))e* ,£ ".- E ».-. 1 ' 

V ^Pi-q 

+ J-^-'K 1+q yb • («0>i + 9>r')7o«(Pi,r)) e"*^"^^ 

V - C/ pi+<? 

V 

(4.56) 

In the "rotating wave approximation" (RWA) the last two terms in Eq. (4.56) 
are neglected, since they describe rapid oscillating processes with frequencies 

uj m 2m/ h. We calculate = — ^[7 , H ion ] using Eqs. (3.32) and (4.56) in the 
RWA- approximation. After some algebra, which is given in Appendix D, we 
can find the expression 
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/„(-*) 



iem r d 3 p i d 3 p f 



I — 

J (2n 



Nfi = v(Pi,Si) tr D 

n { ; ] =u(p f ,s f )ti D 



2h J (2nK) 3 / 2 (2vrft) 3 / 2 \ E p E 

E {K p) dt pS d PvSt + jvw hi h Pit J , 

$f — hty — m 



7o" 



E Pf -hk 



1 ^(9)e* (J VV' 

(4.57) 
^( P/ , S/ ),(4.58) 



^ + ^ — m 



7o 



7 ^ 1 7 ° ~ — 7 M 



E Pf -hk 



Pi +hk 



ufa, Si ). (4.59) 



It is now straightforward to calculate the force-force correlation function, fol- 
lowing the steps in the RPA-approximation (see Appendix C). As seen from 
Eq. (4.55) we need the expression for the real part of the transverse force-force 
correlation. With the definitions Eqs. (3.19) and (3.20) we find after perform- 
ing the contractions using the Wick Theorem [Eq. (C.2)] and integrating as 
well as summing over the variables appearing in the delta-functions, we obtain 



e 2 m 2 _ / ,- j(_ u;+ j J? )t 



-ui+ir, 4h 2 f3 



Re / dte 1 



J di J 



d 3 p i d 3 p f 



(2ixK) 3 (2nhY 



1 (CW)V ( ^^ M 



E P E P , 



E\f(E Pf )[l-f(E p ) 



+ f(E P )[i-f(E P : 



(4.60) 



The integration over t and t can be performed, and we obtain, with the Dirac 
identity lim v ^ +0 (x ± ir])^ 1 = P(l/x) =F i^5(x) 



Re(j x , j A 



m 2 7r r d 3 p i d 3 p f 



-iv AH 2 uj(3 J (2nh) 3 (2itK) 



5(E Pf - E Pi - M 
EpEp. 



«n(^)) 2 (e' 



2/ /3(E P -E P .) 



f -i' 







E /( V [l - KEp)\ + /(E) [l - f(E pi )} (4.61) 

Using the first two relations of Eq. (C.7) as well as Eq. (4.55) and writing 



Ei = E p . and E f = E Pf 



Re(j x , j A 



e 2 m 2 rc r d 3 p t d 3 p f S(E f - - hu) 



r d°p i <n 
J (2nh) 3 (2^ 



-u+i V 4h 2 uj(3 J (2ttH) 3 (2nK) 3 E V E V . 
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e { tew 2 [/w - / w 



+ 



em ( , e) 



/(£/) - /ra 



(4.62) 



Finally with Eq. (4.55) we find for the absorption coefficient 



a(u>) 

E 



lim 



4 2 

emf 



d 3 Pi d 3 p f 5{E f -Ei-huj) 



fc-o 4h 2 u;k 2 J (2ixhf (2nh) 3 



f(E f ) - f(E t ) 



E Pf E p . 



+ 



f(E f ) ~ f{E % 



(4.63) 



It should be emphazised, that the transition matrices Njf^ and Njf^ carry 
a fc— dependence (see Eqs. (4.58) and (4.59)) leading to a finite absorption 
coefficient in Eq. (4.63). As already pointed out the polarization vector ef for 
a linear polarized wave is one of the two transverse modes, labeled here by 
% = 1,2. 

Within the approximation applied here, we observe from Eq. (4.63), that we 
have an electron as well as a positron contribution, responsible for the absorp- 
tion of the external wave. In first order perturbation theory these two terms 
are not coupled to each other. Writing the result in terms of the transition 
matrices Njp and A^ e) , it can further be observed, that our result corresponds 
to the well known Bethe-Heitler [15] formula, if treating the absorption of the 
wave by a single incoming electron with energy Ei and an outgoing electron 
of energy Ef. 

In [16] non-relativistic results for the absorption coefficient in different ap- 
proximations are derived. The expression given in the paper in terms of the 
transition matrix correspond to the non-relativistic result of Eq. (4.63) if the 
matrix element is evaluated in Born approximation. The final result is written 
for the complex collision frequency u, which is related to the absorption coef- 
ficient a by a{uj) = oj^Juj 2 Re v(uj) , with u pl the electron plasma frequency 
and the index of refraction n(uj) = 1. 



5 Discussion and Outlook 



We demonstrated in this work, how the hyperplane formalism can be used 
for a manifest covariant density matrix formulation of relativistic plasmas. 
The one-time description allows to formulate an initial value problem, which 
can lead to considerable simplifications in short-time pump-and probe exper- 
iments. A covariant scheme is developed in the hyperplane formalism, where 
Heisenberg operators are defined on spacelike hyperplanes in Minkowski space. 
In particular, the construction of the quantum Hamiltonian starting from the 



25 



classical QED Lagrangian, making use of the canonical quantization scheme 
is shown. 

From the Liouville von Neumann equation an initial value problem in the 
hyperplane formalism is formulated. The approach used in this work is a gen- 
eralization of Zubarev's method of the relevant statistical operator. For the 
case of moderate fields we applied the linear response approximation, where 
the statistical operator is expanded near its equilibrium solution. From the 
self-consistency relation, which determines the Lagrange multipliers of the 
generalized Gibbsian ensemble, we obtain the response equation, defining the 
susceptibility tensor on the hyperplane in terms of current-current correlation 
functions. The relativistic susceptibility tensor, which displays the response 
of the fermion current to an external electro-magnetic wave, is calculated in 
the RPA approximation. The result agrees with familiar results, published for 
instance by Tsytovich [24]. 

Further it was demonstrated how to include collisions into the formalism in 
a systematic way, making use of perturbation theory. For that reason the 
current-current correlation function is expressed by force-force correlation 
functions by partial integration. The force-force correlation functions contain 
collisions since the force operators are calculated from the von Neumann equa- 
tion with the interaction part in the Hamiltonian. The advantage of the rep- 
resentation of the current- current correlation function, is the appearance of 
the force-force correlation function in denominator, which is convenient for a 
perturbative expansion. 

As an illustration, we derived an expression for the absorption coefficient of 
inverse bremsstrahlung in first order of the force-force correlation function 
(which corresponds to second order in the interaction). We made use of the 
adiabatic approximation, where the dynamic of the positively charged plasma 
component is frozen as well as the Born approximation in calculating the 
correlation functions. The result given here is a generalization of the Bethe- 
Heitler formula for the case the absorption in a electron-positron plasma. The 
interaction is described by the Coulomb interaction. 

Extensions to the approximations assumed here in this work can be done in 
different ways. It is possible to consider higher orders in the expansion of the 
current-current correlation function in terms of the force-force correlation func- 
tion. This will allow to describe higher order processes, like for instance pair 
production. Further interactions can be taken into account, like the electron- 
electron interaction, or the radiation part of the Hamiltonian, which implies 
to couple photons into the plasma. 

A general scheme to derive kinetic equations in the hyperplane formalism is 
explained in [1], valid also for the case of strong external fields. However, 
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processes beyond the RPA approximation are hard to be calculated, since 
coupled kinetic equations for the fermions and photons are to be solved. 



Appendix A 

Commutation relations on hyperplanes 

The constraint equations for the canonical variables A^i and on the hyper- 
plane a n r is written in the form Xn( x ±) = 0> where 

Xl (x ± ) = V^(xJ, x 2 (zj = V M m(xJ, 

(A.l) 

X 3 (x ± ) = n^At(x ± ), Xa(x ± ) = njlt(x ± ). 

For any functionals $ x and <£> 2 of the field variables A ± and U ± , we define the 
Poisson bracket 

a n.r 

where the constraints are ignored in calculating the functional derivatives. 
Applying this formula to the canonical variables we obtain 

[At(x ± ),n ±u (x' ± )] P = S^S 3 (x ± -x' ± ) (A.3) 

with the three-dimensional delta function (2.17). All other Poisson brackets 
for the canonical variables are equal to zero. In the Dirac terminology, func- 
tions (A.l) correspond to second class constraints since the matrix 

C NN ,(x ± ,x' ± ) = [Xn( x ±),Xn>( x '±)]p ( A - 4 ) 

is non-singular. A straightforward calculation of the Poisson brackets shows 
that the non-zero elements of C are 

x ',l 

(A.5) 



C 12 {x ± ,x' ± ) = -C 21 (x ± ,x' ± ) = -V M V^ 3 (x ± - 
C%A x ±i x '±) = ^C^Jx^x'^) = 5 3 (x ± — x' ± ). 
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According to the general quantization scheme [20,21], commutation relations 
for canonical operators are defined by the Dirac brackets for classical canonical 
variables. In our case the Dirac brackets are written as 



[$i,$ 2 ] D = [ $ i>* 2 ] P 

- J da J da' [$, (x ± ,x> ± )[x N ,(0,<5>2] P (A.6) 



'n.r 71, T 



(summation over repeated indices). The inverse matrix, C N 1 N ,(x ± , x' ± ), satisfies 
the equation 

J da C N - N /i(x ± , x ± ) C N /t N /(x ± , x ± ) = 5 NN i 5 (x ± x ± ). (A.7) 



Since the matrix elements (A. 5) of C depend on the difference x ± — x' ± , 
Eq. (A. 7) can be solved for C^ 1 using a Fourier transform on cr , which 
is defined for any function f(x) as 

f(r,pj = J d 4 xe^ 5{x-n- r) f(x). (A.8) 



The inverse transform is 



f(x) = f(r,x ± ) = J ^- 3 e-^5(p.n)f(r,p ± ). (A.9) 



(2k) 



If we perform the Fourier transformation in Eq. (A. 7), we find by insert- 
ing (A. 5) that the non-zero elements of C~ l are 



C£(x ± M = -C^(x ± ,x' ± ) = -j -^e-*M6(p-n)± 

C u \x ± ,x' ± ) = -C^(x ± ,x' ± ) = -5 3 (x ± -x' ± ). (A.10) 

Now the Dirac brackets (A.6) for the canonical variables are easily calculated 
and we obtain 



[A^x ± ),W ± (x' ± )] D = c^(x ± -x' ± ), (A.ll) 
[At(x ± ),A^x' ± )] D = [W ± (x ± ),W ± (x'J] D = 0, (A.12) 

where the functions c^ v (x ± — x' ± ) are given by Eq. (2.25). According to the 
general quantization rules, the commutation relations for canonical operators 
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correspond to i[. . .] D . Thus, in the hyperplane formalism, the commutation 
relations for the operators of EM field are given by (2.23) and (2.24). Obviously 
these relations are valid in the Schrodinger and Heisenberg pictures. 



The anti- commutation relations on hyperplanes 



To find the anticommutation relations for the fermion operators on the hyper- 
plane cr nT , it is sufficient to consider a free Dirac field. Our starting point 
is the standard quantization scheme in the frame where x^ 1 = (t, r) and 

n M = (1,0,0,0) (see, e.g., [18]). In that case the field operators ip a and ip a 
can be written in terms of creation and annihilation operators according to 

J K^) y 2e(p) s =±i 

d 4 p 5{p° - e(p)) 



^ip-x 



jMp) 



E [4(p)««(p)e-*- x +S!(p)«„(p)e*- 



s=±l 



where e(p) = \jp 2 + m 2 is the free fermion dispersion relation. Constructing 

the expression {^ a {x),^ a ,{x')} for two arbitrary space-time points and recall- 
ing the anticommutation relations 

{b s (p),bi,(p')} = {d s (p),dl,(p')} = 5 ss ,5\p-p'), (A.13) 



as well as polarization sums 

l^Pfi + m 

we arrive at 



E U as(P) U a's(p) = 
s=±l 



S V as(PK's(p) 
s=±l 



o?p 1 



(2tt) 3 2e(p) 
+ 



{ [7% + m 



YPf, - m 



Q -ip-{x-x') 
ip-(x-x') \ 

J ' 



where p° = yp 2 + m 2 . Using 
d 3 p 1 



f d"P 1 f d 4 p 2 _ 
J (2tt) 3 2e(p) 7 (2vr)3 



p°>0 



(A.14) 



(A.15) 



Eq. (A.14) can be rewritten in a Lorentz invariant form 
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d A p 



7>„ + ml Q-ip-i?-*') §( p 2 _ 

^ on' v 



m 



p°>0 



+ 



TV a ~ m 



, e^( x - x ') 5(p 2 - m 2 ) 



p°>0 



(A.16) 



The anticommutation relation on the hyperplane a n T is now obtained by set- 
ting x = nr + x ± and x' = nr + x' ± . In calculating the integrals, it is convenient 
to use the decomposition p^ = n^p,, + p±, (p., > 0). Then we get 



d 4 p 5 (p., -e(pj) 



^(r,xJ,4(r,<)} = /^ )3 ^ 

x { biiPy + iIp^ + 



e ^J-m 



+ 



7„P„ + 7>x, - m] e^K-T)} (A.17) 



with the dispersion relation on the hyperplane 

e(P±) = \f-P~Z^- + 



(A.18) 



Finally, changing the variable p ± — > — p ± in the second integral in Eq. (A.17), 
we obtain the anticommutation relation (2.26). The relations (2.27) can be 
derived by the same procedure. 



Appendix B 



The relevant statistical operator in linear response 



In order to rewrite Eq. (3.12), we make use of the operator identity 



and obtain 



H 

\ o 



alz e< c ^ C 2 e~ zC i | 



(B.l) 



^ o 
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X 



/ da E<(^;r)4 At (^)e z/3 ^-^})^oW • (B.2) 



The Lagrange multipliers 0^ are already first order contributions for the de- 
viation from the equilibrium. Therefore in zeroth order, we can express the 
Schrodinger operators Bf(x ± ) by some r-dependent Heisenberg operators ac- 
cording to 



e" 2 ^ B?(x ± , r ) = tf (izf3, r ) B?(x ± , r ) U o (iz0, r ) 

= B$(x ± ,iz(3) , 



(B-3) 



where we assumed, that at r Schrodinger and Heisenberg operators coincide. 
The free evolution operator U is given by 

U q (t, r ; n) = e^ 6 '^'^ . (B.4) 
Making use of Eq. (B.3) in Eq. (B.2) we find the relevant statistical operator 



£ rcl (n,r)= d(7 <t>^{ x ) J dz 

\ in n 



Bg(x ± ,izP) )Q (n) 



(B.5) 



and finally Eq. (3.15) is obtained after Fourier transformation. 



The irrelevant statistical operator in linear response 



In order to express Ag(n, r) [see Eq. (3.4)] in linear response we first consider 
the commutator 



g rcl (n,r),H T (n) = g rel (n,r), H s (n) + H^ xt (n) 



(B.6) 



and note, that the zero order part Q (n), H s (n) = vanishes by definition. 
There are two first order terms 



£o( n ),#cxtH 



and 



(B.7) 



where is the first order contribution of g Tcl . Making use of the Kubo- 
identity 
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C 2 ,e e i] = J dze zd i [d 2 ,C 1 
o 



(B.8) 



we find 



Since 



3 U ,Q 



n 



(3 J daA^ xt (x) J dze-^-^) 

°n 

X,(H s -vQ)]e z ^-^) Qo (n). (B.9) 
H s: Q = (Q is an integral of 



3 [ii J a 



= and also 



motion) we can write 



g (n),HUn)} = P J daA^(x) J dze~ z ^ [j^H,] e z ^g (n) . (B.10) 



Since the external field A^ xb (x) appears in (B.10) the commutator j^,H s 
only treated in zero order, and we can write 



is 



9 T Jf, = -' 



Jui H — —i H s 



Making use of Eq. (B.3), we find 



i 



0o( n )iHLt( n )] =i (3 J daA% xt (x) J dzj^x ± ,i(3z 



i > I <i(TA.^ t {x) I az]^x^ipz)Q [n) 

°n 
1 

ipe~ lk w T J dz AZ xt (k)j^(k ± ,if3z)e (n) , 
o 



with the last expression written in Fourier space. 



(B.11) 



(B.12) 
(B.13) 



The second commutator in Eq. (B.7) can be written in linear response by 
making use of Eq. (B.8) 



Q^(n,r),H s (n)\ = i/3£ / da ^(x) I dzB^x ± ,zzP)g (n) . (B.14) 



In Fourier representation Eq. (B.14) reads 
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g^(n,r),H s (n)\ = i/35>£(A;) / dz B^(k ± ,izP)g 



n 



(B.15) 



The derivative term in Eq. (3.4) is calculated using Eq. (3.15) 



_L 

d T g Icl (n,r) = ~t(3k l} e- tk f T f dz B^(k ± ,iz/3) g (n) . (B.16) 



Constructing the expression for Ag(n, r) in the linear approximation [Eq. (3.4)], 
the evolution operators U are to be taken in zeroth order and yield a r- 
translation. From the Eqns. (3.4), (B.13), (B.14), (B.16) one finds 



Ag(n,T) = -P J dr'e- r ' (T - T '^- ik \\ T ' J dz 

-oo 

]T [B^(r' - r + iz(3) - ik^Yir' - r + 

(B.17) 



Eq. (B.17) can be rewritten using the transformation r = r — r' and Eq. (3.16) 
is obtained. 



Appendix C 



The RPA susceptibility tensor 



We evaluate the current-force correlation function in the RPA approximation 
with Eq. (3.32), (3.34) as well as j v {k) = j u (—k) and the short-hand notation 

i Jvify) u)+ir) = U/j, 'ji>)ui+ir) 



k 



oo & 

l„ ■ lu}^ = 'die^ + ^ J 'dtTr{j\(k,t-iht)l(-k)g } 



o 







• 2 2 

ie m 

hp 



J dt e^ +ir ^ J dt J d 3 p d 3 p' 



- E 
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I — ' I 1 1 

Trj •d ps d p _ lrnks ,-.:d p , r d p ,_ hkr ,: g \ 

/-fife 



tr D {(v(p + hk, s'^vfa s)) (u(p' - /&),r')7„u(p',7-))} 
*y - ^-fife] 



+ 



1 



1 



I 1 



tr D { (u(p - hk, s')7^(p, s)) (u(p' + hk, r')-f v u(p', r)) } 

\e,-e, Jp-l^-v^x*-^ 



+ 



1 ~ L ~ \ 

\ ,- \ p Trj :d_ p+nks ,b ps ::d p , r b_ p , +hkr t: g > 

tr D { + s')7„«(p, s)) + ^A;, r'favtf , r)) } 



+ 



e -i(^- s P - M )(*-^") 
r 



I , C L 

Tr| :^l p _aA:,s'^P,s : : d-p>-hk,r'K',r : £o} 



\ -Ep+fife \ Ep'+hk 
tr D { - ft/c, s')t^(p, s )) - r'favtf, r)) } 



p'+fife 



(C.l) 



The notation tr D {. . .} stands for the Dirac trace. Making use of Wick's theo- 
rem [18,26], the non- vanishing contractions, which are denoted by the brackets 

b ps b p , s ,, can be written in the form 



p')S ss '(l-f(E p j) 

p')t ss >{l-f(E p )) ■ (C.2) 

Further we introduced the Fermi-distribution functions / and / for particles 
and antiparticles respectively 

f( E p) = e P{E p -») + l > f( E p) = e P{E p+ ») + I ( C - 3 ) 

It should be noted that contractions inside the normal order do not contribute 
and that crossing contraction lines (the third term in Eq. (C.l)) give an extra 



Tr{&J s ,& pV £o} = Kp- 

Tr {v^i'^o} = s (p- 

Tr {dl s ,d p , s ,Q } = 5{p- 
Tr {d ps ,dl, s ,g } = 5(p - 
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minus sign. Eq. (C.l) can now be written in terms of Fermi functions if the 
integration over p' and the summation over r and r' is performed 



LO+irj 



■ 2 2 

%e m 

hp 



J J J {2nhf E p 



^p+hk 



p+hkj 



M E v- E v+h k )^- ihl ) 



tr D { (v(p + hk, s')7 M u(p, s )) ( v (p, s)lyv(p + hk, s') } 



+ 



E, 



p—hk 



f(E p _ hk )\l-f(E p 



Ep-hk ~ E p 



tr D { (u(p - Hk, s')7„«(p, s )) ( M (P> «hXP - hk, «')) } 



p—hk 



l~f(E p _ hk ) l~f(E p ) 



E p + E p _ hk 



+ 



ti D {(y(-p + hk, s')7 M «(p, s)) (u(p,s)7 I/ u(-p + ^A;, s'))} 
p f(E p+hk )f(E p ) 



J p+hk 



E p + E p+m 



M E v+ E P+hk )$- if > i ) 



tr D { (u(-p - hk, s')~f^v(p, s)) (v(p, s)^ v u(-p - hk, s')) } [ . (C.4) 



Finally the spin summation rules 



J2u(p,s)(g)u(p, s) = —(p' + m) 
Y, V (P> s)®v(p,s) = ^(P 1 -™) 



(C.5) 
(C.6) 



and the integrations over t and i can be performed in Eq. (C.4). With the 
relations 



f(E p )(l - f(E p+hk )) (e+^-W - l) = -(f(Ep) - f{E p+m )) 

f{E p _ hk )(l - f(E p )) (e-^-V«) - l) = /(£y - J 

(l " f(Ep)) (l - /( (e-^-W - l) = -1 + + /(Vtt) 

W( e+/,(Bp_W - X ) = 1 - W " /(W ( C - 7 ) 
the equation (3.35) is obtained. 
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Appendix D 



Calculation of j^(k) 



From the Eqs. (3.32) and (4.56) we find 



/„(-*) = 



d 3 q l3 l3 lm I m 
drpdrpxJ—.l—A? 



m 



h J (2nhf 

•<7 /ft • -H (ft • 

• u p-hk,s' u p,s' 1 • u p 1 -q,r' u p 1 ,r' 



EVE 

P V Pi ss'rr' 



Xq) E { 



m 



l( E p- E p-hk +E P 1 ~ E P 1 -g) t 



+ 



m 



m 



+ 



m 



m 



y Ep—fik y Ep^—q 

x (u(p + s')7 M u(p, s )) (u(Pi - 9, 07o u (Pu r )) 

'.ft h ■ -ft h •] p^~ E P +E r+™~ E *>i +E P 1 +J t 

•°p+hk,s'°p,s- > •°p 1 +g,r'°p 1 ,r*J e 

x (u(p + s'^uip, s)) (u(p! + g, r')7o w (Pu r)) 

V + ^ fc y _E |_ p -nfc ' s P> s " 5 • a p 1 -q,r' a p v r-\ « 

x - t&, s'h„u(p, s)) (u(p! - g, r')7o u (Pu r )) 

£ . .£t £ .1 p i^ E P^ E p+hk~ E V 1 +E p +q) 1 

• u -p-hk,s' u p,s' j • u p 1 +g,r' u p 1 ,r*J c 

- Hk, s')t m «(p, s )) («(Pi + 9, r ')7o M (Pi> r )) 

y -L^p—tik V I - J p 1 —q 

x + s')7 M u(p, s )) (u(Pi - 9, 07o u (Pu r )) 

e S(- E P+- E p-fifc - - E P 1 +- E P 1 +9)* 



Ep+fOt y E Pi+q 



x 



+ 



m 



Ep—hk y Ep^+q 



•ft (ft • -ft h 

• u -p+hk,s' u p,s' J •°'p 1 +(j,r /U 'p l ,r* 



X 



(u(-p + s')tXp, s )) + 9> 07o u (Pu r )) } • 



(D.l) 



After calculating the commutators in Eq. (D.l), the integration over p 1 as well 
the summation over the primed variable can be performed. Finally, using the 
spin summation Eq. (C.5) and (C.6), we find the following equation for the 
force operator (only the non-vanishing terms in the RWA- approximation are 
kept) 



'■i i \ «em f d 3 q 3 „ 
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E r V E 



e Ji( E p- E p 



-q — hk ' 



p—q—Hk 



X 



{y{p -q-hk, s^To - ^ — —%v(p, s)) d ] ps d p _ q _ hksl 

^p—Hk 



A 



i E v+ a - E v -hk) t 



E 



p+q V Ep-hk 



Qh\^p+q 



X (u(p - fife, S')7 M — j^loV(P + Q, sj) dl+qJp-Hk,s> 



+1 



Ep-q 



Qh( E p+hk 



E. 



p+hk 



x(u(p + hk, s'h M — ^7o^(P - 9, s )) ^+wfc, s ' V, )S 




-Ep V Ep+q+hk 



e -Ji( E P~ E p+q + hk) t 



X 



[u(p + q + hk,s )7 7 u(p, s) J ^ 

&p+hk 



J p+q+hk,s'^ > p,s 



(D.2) 



We perform the shift p — > p + hk in the first term and p ^ p — hk in the last 
term of Eq. (D.2) and introduce the initial and final momentum variables, p i 
and according to 



Pi = p — hk , Pf = p + q 
Q = Pf-Pi-frk , " ' ' . (D.3) 

Pi=P-Q , p f =p + hk 

Using these redefinitions we can obtain Eqs. (4.57) - (4.59) in terms of the 
transition matrices Njf^ and Nffi for electrons and positrons respectively. 
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